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Abstract 



Small size QCD instantons may spoil the axion solution to the strong CP 
problem if QCD is not asymptotically free at high energy scales. We exam- 
ine this issue in supersymmetric models using a manifestly supersymmetric 
scheme to compute the axion potential induced by small size instantons. Ap- 
plying this scheme for a class of illustrative models, it is found that the result- 
ing high energy axion potential is highly model-dependent, but suppressed 
by more powers of the soft supersymmetry breaking parameters and/or of 
the other small mass scales than what is expected based on a naive instanton 
graph analysis. Our analysis suggests that the axion solution is stable against 
the small QCD instanton effects in a wide class of supersymmetric models 
even when QCD is not asymtotically free at high energy scales. 
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1. INTRODUCTION 



One of the most attractive solutions to the strong CP problem is to introduce a 
spontaneously broken U(1)pq symmetry whose explicit breaking is entirely given by the 
QCD anomaly 0. Once such a global U(l) pq is assumed, it is usually taken for granted that 
@qcd is dynamically relaxed down to a sufficiently small value by the vacuum expectation 
value (VEV) of the associated pseudo-Goldstone boson, the axion. However there are some 
class of models that this may not be true. If QCD is not asymptotically free at energy scales 
above the weak scale, small size QCD instantons may induce a significant high energy axion 
potential Vhe whose minimum does not coincide with the minimum of the conventional low 
energy axion potential Vle- If Vhe ~ 10~ 9 Vle, the dynamical relaxation mechanism of Qqcd 
would be spoiled and then one looses the original motivation to introduce U(1)pq. 

Perhaps the most interesting class of models in which the effects of small size instantons 
can be sizable are supersymmetric models. Some supersymmetric Grand Unified Theories 
(GUTs) like Eq and also many supersymmetric models derived from superstring theories 
contain extra colored multiplets which would dramatically change the QCD f3 function at 
energy scales above the weak scale, even making it positive in some cases. In fact, small size 
instanton contribution to the axion potential in models with extra colored particles has been 
studied before based on a simple dimensional analysis for instanton graphs. However, 
when applied for supersymmetric models, this scheme can overlook a possible cancellation 
between different instanton graphs which is due to supersymmetry (SUSY), thereby yielding 
a highly overestimated axion potential. To avoid this difficulty, in this paper we introduce a 
manifestly supersymmetric scheme to compute the axion potential and apply it for a class of 
illustrative models. Although an explicit analysis is made only for a limited class of models 
and also the resulting axion potential is highly model-dependent, our analysis suggests that 
the effects of small size instantons in supersymmetric models are much more suppressed than 
what is expected based on a naive instanton graph analysis. This implies that the axion 
solution to the strong CP problem is stable against the effects of small instantons in a wide 
class of SUSY models even when QCD is not asymtotically free at high energy scales. 

To proceed, let us discuss in more detail how the dynamical relaxation of Qqcd can be 
spoiled by the high energy axion potential induced by small size QCD instantons. At the 
axion scale f a where U(1)pq is spontaneously broken, the axion effective lagrangian is given 
by 

£axion = \{d,af + ifyaJ" + (1) 

2 fa 62TT Z f a 

where denotes a generic current in the model, and F and F are the gluon field strength 
and its dual, respectively. Integrating out the degrees of freedom at energy scales above the 
weak scale Myy, the nonderivative axion coupling to the QCD anomaly, i.e. aFF, would 
induce a high energy axion potential Vhe through the QCD instantons with size p <C M^ 1 : 

VkE = -V cos(a// a - 7 ). (2) 

Both the size of this axion potential (Vo), and the location of the minimum (7), depend upon 
the details of the physics at high energy scales around p~ x . Unlike Vhe, the low energy axion 
potential induced by the dynamics below the weak scale can be unambiguously computed. 
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Applying the chiral perturbation theory for the low energy QCD instanton amplitudes, one 
finds 

Vle = -flml\Jm 2 u + m 2 d + 2m u m d cos(a/ f a )/(m u + m d ), (3) 

where m u and denote the current masses of the up and down quarks, f n = 94 MeV and 
m v are the pion decay constant and mass, respectively. Here the axion VEV ((a/ f a )) is 
defined to correspond to the QCD vacuum angle which is required to be less than 10~ 9 from 
the non-observation of the neutron electric dipole moment: 

e QCD = (a/ f a ) * 10- 9 . (4) 

The axion VEV is then determined by the total axion potential V ax ion = Vle + Vhe- Obvi- 
ously, in order for the axion VEV to satisfy the above phenomenological bound, the high 
energy axion potential has to be suppressed as V^j f^m\ ~ 10 _9 7 _1 . If one wishes to achieve 
CP-conserving QCD independently of the details of CP violation at high energy scales, i.e. 
independently of 7, one would need 

x = Vhe < , /_v 

In fact, there is a negligibly small extra low energy potential || 

5V LE ^ 10 _14 /^m^sin^ M sin(a// a ) 

due to the low energy dynamics involving the CP violation by the Kobayashi-Maskawa phase 
8 km- This would shift the axion VEV also, but the resulting shift (Sa/ f a ) — 10~ 14 sin 5xm 
is well within the experimental bound. 

In the next section, we describe how the instanton-induced Vhe can be computed in 
generic supersymmetric models in a manifestly supersymmetric way. In sect. 3, we apply 
this method for a class of illustrative SUSY models including extra quark multiplets and 
make a numerical estimate of the largest possible Vhe- Based on this analysis, we conclude 
in sect. 4 that the axion solution to the strong CP problem is stable against the effects of 
small instantons in a wide class of SUSY models even when QCD is not asymtotically free 
at high energy scales. 



2. SMALL INSTANTON INDUCED AXION POTENTIAL IN SUSY MODELS 

In this section, we discuss how the high energy axion potential induced by small size QCD 
instantons can be computed in a manifestly supersymmetric way. Generic supersymmetric 
models contain a SUSY breaking sector which would provide a dynamical seed for SUSY 
breaking. The resulting SUSY breaking is transmitted to the observable sector through 
some messenger interactions characterized by the messenger scale M m . Integrating out the 
fields in SUSY breaking sector and also in the messenger sector, one is left with an effective 
lagrangian of the observable sector fields in which SUSY appears to be softly broken. This 
effective lagrangian describes the observable sector dynamics at energy scales below M m . 
The size of the messenger scale M m depends upon the character of messenger interactions. 
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In gravity- mediated models 0, M m is given by the Planck scale Mp, while in gauge- mediated 
models |7j M m is somewhat model-dependent but typically far below M P . At any rate, soft 
SUSY breaking in the observable sector is operative only at energy scales below M m . As 
a result, QCD instantons with size p ~ M' 1 do not contribute to the axion potential due 
to SUSY fermion zero modes, e.g. gluino zero modes. To be definite, here we concentrate 
on gravity-mediated models, however our study can be easily applied for gauge-mediated 
models also. 

Our starting point is the effective lagrangian of observable sector fields at energy scales 
around p -1 where p denotes the size of the small QCD instantons: 

C = J dWOrj^Mi + ^iJ d 2 6YW a W a + h.c.} + [J d 2 6W(<$> h Z) + h.c]. (6) 

Here $j denote generic chiral superfields in the model, W a are the spinorial chiral super- 
fields for the gauge multiplets, rj is a real spurion superfield whose auxiliary .D-component 
corresponds to the soft scalar mass, and finally Y and Z stand for the chiral spurions whose 
auxiliary F-components correspond to the gaugino mass and the soft A (or B) coefficients, 
respectively. Explicitly, we have 

r] = i + m 2 e 2 e 2 , 

8tt 2 „ „ 

Y = — + iO QC D + 167i 2 m 1/2 9 2 , 

Z = {Z x = (l + A6 2 )\, Z^ = (l + B9 2 )p}, (7) 

where A and p denote the trilinear Yukawa couplings and the bilinear //-parameter in the 
superpotential. It is assumed that all soft breaking terms are described by the universal soft 
scalar mass mo, the universal gaugino mass mi/2, and the universal trilinear A and bilinear 
B coefficients, which are presumed to have the weak scale values. The size of the small 
instanton-induced axion potential is not sensitive to the detailed form of the soft parameters 
and thus its order of magnitude can be estimated in the scheme assuming a universal form 
of soft parameters. 

Once the effective lagrangian at energy scales around p _1 is given as @, the low energy 
interactions induced by small QCD instantons can be summarized by [H 

£ ins = J d 2 6d 2 6 [e^ Y K cS (<S> u $*; D a , D a , Z, Z*; 77) + h.c] 

+ [ J d 2 6e- n * Y W cS {^ Z) + h.c], (8) 

where n\ and n 2 are positive integers corresponding to the instanton winding number and 
D a denotes the supercovariant derivative acting on either superfields or spurions. In order 
to generate an axion potential, all external fields with vanishing VEVs, e.g. external quark 
superfields, in the instanton amplitudes have to be integrated out to induce a term including 
only the superfields with nonvanishing VEVs. Then the axion potential can be derived from 
the terms in W e s or K e g by replacing all the superfields by their VEVs, the spurion Y by 
Y + ia/f a , and finally integrating out the Grassmann coordinates 6 and 9. Note that the 
supercovariant derivatives are not allowed to appear in W e s since Z)$ = DZ = and also 
D 2 Q {D 2 Z) is not a chiral superfield (spurion). 
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Let m so f t collectively denote the soft mass parameters, i.e. 

{^soft} = {m ,m 1/2 ,A,B}. 

Usually any nonzero VEV of the auxiliary F-component of the observable sector field $j is 
due to the presence of soft SUSY breaking terms and thus it is suppressed by some powers 
(typically a single power) of m soft , e.g. 

(^^(l + rrW? 2 )^), (9) 

where <pi is the scalar component of the superfield $j. One simple consequence of this 
observation is that the axion potential arising from K e s is suppressed at least by two powers 
of pm so f t since it is obtained after the Grassmann integration over d 2 8d 2 8, while the axion 
potential from W e g is suppressed at least by a single power of pm so ft- 

The instanton-induced effective Kahler potential K c r and the effective superpotential 
Weg can be systematically computed by expanding them in powers of the small dimensionless 
quantities: p$j, pD a D a , Z\, and pZ^. In fact, such an expansion is severely constrained 
by the selection rules dictated by the internal symmetries of the tree level superpotential 
W($>i, Z) in (^|) under which both the superfields and spurions transform. Another useful 
requirement is that it should have a sensible limiting behavior when some of the couplings 
are turned off. As we will see, together with the holomorphy, these requirements forbid W e g 
to include a term which would contribute to the axion potential by replacing the superfields 
by their VEVs. As a result, the leading contribution to the axion potential comes from the 
effective Kahler potential K e $ and thus is suppressed at least by two powers o/pm so f t . 

Once the superfields in K e $ are replaced by their VEVs and also the Grassmann inte- 
gration is performed, the expansion in powers of pD a D a , pZ^, and Z x appear to be an 
expansion in (pfa), pm so f t , pp and also the Yukawa couplings in the model. Generically the 
field VEVs are calculable in terms of m so ft and/or the other mass parameters, e.g. p and 
Mp, in the model. One can then easily identify the dominant term by counting the powers 
of pm so f t and pp (and also of the Yukawa couplings if necessary) in the resulting axion po- 
tential. Although the selection rules determine how the high energy axion potential depends 
upon the various parameters and the VEVs, they do not fix the additional suppression fac- 
tor arising from the loops in the instanton graph. The axion potential induced by small 
instantons is essentially given by closing the fermionic zero modes in the instanton-induced 
multiple fermion operator in the dilute gas approximation. Once the dominant term in the 
expansion of K e $ were identified, the corresponding instanton graph can be identified also. 
This allows us to estimate the loop suppression factor by counting the number of loops in 
the graph, thereby completing the estimate of the size of the high energy axion potential. 

The above procedure computing the instanton-induced axion potential is based on a man- 
ifestly supersymmetric formulation and thus automatically takes into account the possible 
cancellation dictated by SUSY. Anticipating SUSY cancellation between different instan- 
ton graphs would be highly nontrivial in a naive diagrammatic analysis. In fact, in many 
cases one can draw an instanton graph which is consistent with the selection rules and the 
sensible limiting behavior, but with a less powers of pm so f t than those required by SUSY. 
This implies that such graph is cancelled by other instanton graphs which are related to the 
original graph by SUSY. Our method automatically avoids such misleading graphs since it 
includes only the instanton graphs allowed by both SUSY and the selection rules. Following 
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the recipe described above, in the next section we will compute the small instanton-induced 
axion potential in a class of illustrative SUSY models. 



3. COMPUTATION AND NUMERICAL ANALYSIS 

There are many possible supersymmetric extensions of the Standard Model, starting 
from the minimal supersymmetric standard model (MSSM) to other models including ad- 
ditional matter multiplets. Here we are mainly interested in models with additional colored 
multiplets which would change the QCD (3 function at high energy scales, thereby leading 
to a sizable high energy axion potential generated by small QCD instantons. To examine 
this possibility, in this section we consider a model containing A^-flavors of isosinglet quark 
and antiquark multiplets (h + h c ) and also a gauge-singlet superfield (S) in addition to the 
3 generations of the MSSM quark multiplets (isodoublet quarks Q and isosinglet antiquarks 
jjc _|_ jjc^ This model is a rather simple generalization of the MSSM, but still shows many 
of the essential features of the axion potential induced by small instantons in generic SUSY 
models. 

The tree level superpotential of our model includes 

Ok+l On+1 cl+3 

W = Z U QH U U C + Z d QH d D c + Z h — T hh c + Z 1 -—H u H d + Z 2 — T , (10) 

Mp Mp Mp 

where /, k, n are non-negative integers, and H u and H d denote the isodoublet Higgs super- 
fields. Terms involving lepton superfields are irrelevant for the analysis of the axion potential 
and thus are ignored. Here Z Uyd are 3x3 matrices, Zh is Nh x Nh matrix, and all couplings 
are spurions which are determined by the Yukawa couplings and the universal ^-coefficient: 

Z u = X u (l + A6 2 ) , Z d = X d (l+A9 2 ) 

Zi = Ai(l + AO 2 ) , Z 2 = \ 2 (l + M 2 ) , (11) 
Z h = X h (l+A6 2 ). (12) 

A characteristic feature of the models considered here is that all mass scales are determined 
in terms of the two basic mass scales m so f t and Mp and also the three non-negative integers 
(l,k,n). (See Eqs. (14) and (15).) It is straightforward to repeat our following analysis for 
other type of models including a mass scale which is not related with m soft and M P . 

To proceed, let us first consider the ground state configuration of the model. We assume 
that S has a sizable renormalizable Yukawa coupling, for instance the Yukawa coupling 
Xh with h + h c in the case of k = 0, or a renormalizable Yukawa coupling with lepton-like 
superfields in other cases with k > 0. This would make the running soft mass of S to become 
negative at energy scales around (S). Then the effective potential of S is schematically given 
by 

o2(/+2) 

V « ~m 2 o{t S 2 + -j^g-, (13) 
where we have omitted the coefficients of order unity. The resulting VEV of S is 
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(s)~M P (^y\ (14) 

yielding the supersymmetric masses and fih) of H u + H d and h + h c as follows: 

. (S) n+1 fm soit \7ri 

A* = A *^T - ™** [-j^) , (15) 

where we have assumed that Ai,A2 and A^ are of order unity. These masses can not be 
significantly smaller than m so f t which is presumed to be of order the weak scale, and thus 
the non-negative integers I, k,n are required to satisfy 

n<l, k<l. (16) 

For the case with n — l, the Higgsino mass fin is of order m so ft, and then H u and H d can 
be interpreted as the MSSM Higgs doublets which have the weak scale VEVs generating the 
masses of W , Z, and also of the quarks and leptons: 

(H u ) « (H d ) « m soft for n = I. (17) 

In this case, the observed quark mass spectrum gives 

det(A u A d ) w 10" 16 for n = l. (18) 

However for I > n, we have \xn 3> m S oft an d then if u and should be interpreted as 
additional massive Higgs doublets other than the MSSM Higgs doublets. Such additional 
(superheavy) Higgs doublets appear quite often in SUSY E 6 models and also string orbifold 
models |J. At any rate, for p, H ^> m soft , H u and H d have vanishing VEVs: 

(H u } = (H d } = for l>n, (19) 

and then X u and A^ are not constrained at all. As a final remark on the vacuum configuration, 
one can explicitly confirm that 

I dW\ 

= )« merit (20) 

for the models under consideration, which would assure the expression ([|) for the superfield 
VEVs. 

Given the above described vacuum configuration, the flavors of additional isosinglet 
quarks (h + h c ) affect the QCD (3 function at energy scales above particularly make it 
positive if Nh > 3. Obviously we are then interested in small instantons with p ~ /i^ 1 . To 
be definite, we will consider here only the three distinctive values of A^: the case with one 
h + h c in each generation (Nh = 3), the case with two h + h c in each generation (Nh = 6), 
and the case with three h + h c in each generation (Nh = 9). We then require that the QCD 
fine structure constant a does not blow up at energy scales below Mqut = 2 x 10 16 GeV for 
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Nh = 6 and Nh = 9 which is a conservative condition. We do not require no blowing up of 
a above M GUT to M P because there is a possibility that enlarged gauge group at M GUT can 
change the running behavior of a. This leads to the following conditions on fi^ and thus on 
k and 

N h = 6 : fi h > 1(T 5 (m 8aft M GI7T ) 1 / 2 -> < ^, 

iV h = 9 : fA h > 5 x 10- 4 (m soft M^ r )V3 _ *±I < g. (21) 

Within the range / < 2 which we will concentrate on, all possible sets of (I, k) are restricted 
to (1, 0),(2, 0),(2, 1) for N h = 6 and (1, 0),(2, 0) for N h = 9 by the above constraints. 

As we will see, the axion potential induced by small instantons is highly sensitive to the 
Yukawa couplings Aj (i = u, d, 1, 2, h). For Aj, A2, and A^ which are unknown in any case, we 
simply assume that they are all of order unity. Note that this corresponds to a conservative 
choice when applied to check whether the resulting high energy axion potential satisfies the 
strong CP condition (Q). About X u and A^, we have two possibilities. For the case with 
n = I, H u and H d are the usual MSSM Higgs doublets and thus det(A u Ad) ~ 1CT 16 . However 
for I > n, H u and H d are additional (superheavy) Higgs doublets with vanishing VEVs, and 
then A u and A^ can be assumed to be of order unity once again as a conservative choice. 

Small instanton-induced axion potential has the tunnelling factor e for the single 
instanton contribution (ni = n 2 = 1) which gives the most dominant effect. This tunnelling 
factor can be determined by the one-loop (3 function: 

m = E*L = -U.E)jL, (22, 

where b (E) = 3 for E < fi h , but b (E) = 3 — N h for E ^ We then have 

exp(-2vT/a(p)) w (m soit /fi h ) 3 (p/j, h ) 3 ~ h exp(-2n/a(m soit )) 

w 2 x 10- 31 (pm soft ) 3 ~^(^soft/Mp)^r^, (23) 



where we have used the expression (15) for fi^ together with the numerical value of 
exp(— 2tt / a(m soit )) s;2x 10 -31 for m soft w 1 TeV. 

Another important ingredients in computing the instanton-induced axion potential are 
the selection rules dictated by the underlying theory at high energy scales around p -1 . ^From 
the tree level superpotential, we can find various global symmetries under which both the 
superfields and the spurions transform. Low energy effective theory obtained by integrating 
out high momentum modes should preserve these symmetries. The global symmetries of the 
tree level superpotential ( |10D include 

G = SU(3) Q x SU(3) UC x SU(3) D c x SU(N h ) h x SU(N h ) h c 

xU(l) A x U(l) x x U(l) x > x U(1) R (24) 

where SU(N)^. denotes the SU(N) rotation of the iV-flavors of $j, and the £7(l)-charges of 
all superfields and spurions are depicted in Table 1. 
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Let us consider a G-invariant effective superpotential e~ Y W e s with a sensible limiting 
behavior, which would lead to an axion potential once the superfields are replaced by their 
VEVs. For a small instanton size p < (S)^ 1 , it takes the following form: 



~ Y W cS ~ p- 3 e- Y det(Z u Z d Z h )Z^(p 2 H u H d ) N *(pS) N *Z^ , (25) 



where 



(pMp)*' 1 (pM P )"' ' (pM P ) 



The integers iVj (i = 1, 2, 3, 4) in the above are required to be non-negative to have a sensible 
limiting behavior. This condition can not be compatible with the [/(l)p-selection rule 

^ + ^ + 2 = 0, (27) 

and thus no such an effective superpotential is allowed. 

Since it does not have to be holomorphic, a generic instanton-induced Kahler potential 
which would lead to an axion potential can be written as 



- v , . 7 7 JZA* 1 ( P 2 H u H d \ N2 ( P S\ NS (ZA N4 (pD a D"\ N * 
e K eS = p e aet(Z u Z d Z h ){ M ) ^ j ^_ ] , 



(28) 



where -F(^) is a function of the real spurion r\ = l+ml9 2 9 2 and the supercovariant derivatives 
are understood to be applied for either the superfields or the spurions. The G-selection rules 
require first of all the factor det(Z u Z d Z h ). One may then insert either Ni -powers of Z 1 or 
its complex conjugate, iV 2 -powers of H u H d or its complex conjugate, and so on, to make 
K e g to be fully G-invariant. To make the notation simpler, in the following we will use 
the notation in which a negative power corresponds to the positive power of the complex 
conjugated fields or spurions. In this notation, the above Kahler potential can be written as 

e~ Y K cS = p- 2 e- Y det(Z M ^^)ZfHp 2 ^^) iV2 (P^) JV3 ^ 4 (P^ 2 ) JV5 ^(^), (29) 

where now {Ni} can be negative integers but still are constrained by the selection rules of 
U{l) x , U(l)x> and U(1) R as follows: 

N 2 = N 1 + 3, 

N 3 = (n + 1) iVi + + 3)iV 4 + (k + l)N h , (30) 
N 5 = N 1 + N 4 + 3. 

In fact, there can be two additional suppression factors that are encoded in F{rj). The 
first one is the loop suppression factor which is always there, while the second one appears 
only when p < f~ l . Let 0pq denotes the scalar field responsible for the spontaneous U (1)pq- 
breaking. We then have 

(<Ppq) = (31) 
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where the positive integer Nqw corresponds to the axion domain wall number. Then the 
L r (l)pQ-selection rule, i.e. the invariance under 



9qcd — ► Oqcd — aN DW , (32) 

requires that the axion potential induced by a very small instanton with p < f~ l includes 
the factor (p((I)pq)) Ndw ■ This factor is replaced by the simple phase factor e %a ^ a for larger 
instantons with p > f~ l . Summarizing these, we have 

Ffo = 1) « e. (JL^ (33) 



where 



, ^i(pf a ) NDW for p<f~\ 

€a ~\l for P>f~ l (34) 

and L denotes the number of loops in the instanton graph. 

After replacing the superfields in K e R by their VEVs and also integrating over the Grass- 
mann variable, we find 

Vhe(p) « p~ 4 F(l)e-^(pm soft ) 2+ l iV5 l(pMp)- iVo (p 2 (i/ u ^))l JV2 l 

x (p(S))^det(X u X d X h )X[ Nll x[ Nil , (35) 

where 

N = kN h + n\Ni\ + l\N 4 \. (36) 

Note that the above axion potential includes the model-independent SUSY suppression factor 
(pm so f t ) 2 in addition to the tunnelling factor e~~> and also the model-dependent suppression 
factors due to small field VEVs and couplings. Using the field VEVs fll4|), and (1P9|) 
together with (|2"3] ) and (P3"|), we finally arrive at 

Uhe(p) « 2 x 10- 31 ea p- 4 (pm soft ) 5+ ^-^(pMp)- JVo+ ml 7V3 l 

x (^J (^-) i+1 det(A u A d A ft )A^ l! A 2 (37) 

where 



2|iV 2 | + ^iVs! + |iV 5 | for Z = n, 
iV 5 = { (38) 
^|AT 3 | + |iV 5 | for Z>n. 

Here the loop number L can be determined by looking at an explicit instanton graph yielding 
the above form of axion potential. Note that (H u ) ~ (Hj) ~ rn so { t only for I = n. For other 
cases with I > n, we have = (-Hd) = and thus only N 2 = can yields a nonvanishing 
axion potential. 
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Given the values of Nh = 3,6,9 and (n,l,k) satisfying the constraints (|16D and (pT]) . 
we can now pick out the values of {iVj} (i = 1 ~ 5) which would give the most dominant 
contribution to Vhe for a fixed value of p. It turns out that for all cases studied here such 
values of {N} lead to a negative total power of p in Vhe and thus smaller instantons give a 
larger contribution. Since we wish to see whether Vhe is small enough to satisfy the strong 
CP condition (H), we are interested in the largest possible value of Vhe, i-e. an upper bound. 
We thus take p ~ 1/Mp, being the smallest possible value of p, and also did not take into 
account the / a -dependent suppression factor e a ~ (pf a ) NDW in our numerical study. It is 
then clear from eq. ( |3"7| ) that, for the fixed values of (N h , n, I, k) and pM P « 1, the largest 
possible value of Vhe is obtained for the minimal values of N$ and also of the loop number 
L. In each case, it is rather straightforward to find the minimal value of iVg and also the 
corresponding instanton graph with a minimal number of loops. (See Figures 2, 4, 5, 6.) The 
results are summarized in Table 2 for all cases with I < 2. It turns out that Ng < 2 in all 
cases with I = n <2 and thus N 2 = even when / = n. We stress that the values of Vhe i n 
Table 2 are obtained for the smallest possible value of p (~ Mp 1 ) and the conservative choice 
of the loop factor ^), and also the possible additional suppression factor e a « (pfa) Now 
for f a < was not taken into account. In this sense, they could be largely overestimated 
and thus should be understood as a kind of upperbound on Vhe- 

As we have pointed out, a naive instanton graph analysis can yield a misleading result 
with less power of pm so f t than what is required by SUSY and the internal symmetry selection 
rules. Such a graph should be cancelled by other instanton graphs which are related to 
each other by SUSY. Avoiding this complication was the main motivation to introduce a 
manifestly supersymmetric scheme to compute the small instanton-induced axion potential. 
To make the motivation more clear, let us consider some examples of misleading instanton 
amplitudes. For the case of (n, /, k) = (0, 1, 0) and N h = 9, there can be a diagram like Fig. 
1 which is consistent with the selection rules of all internal symmetries, but includes only 
one insertion of pm so f t (the A parameter in this case). The size of the resulting amplitude 
is estimated to be (for p _1 « / a « Mp) 



which would be much larger than 10 _9 /^m^ if all the Yukawa couplings are taken to be of 
order unity. (Note that I > n in this case and thus H u and H d should be interpreted as the 
superheavy Higgs doublets which have a vanishing VEV and thus can have a large Yukawa 
coupling with the standard model quarks.) Based on this result, one may conclude that 
the axion solution to the strong CP problem can be spoiled by small size QCD instantons 
in this model. However to be consistent with both the SUSY and the internal symmetry 
selection rules, the instanton amplitude should include at least two powers of pm solt and thus 
the above amplitude should be cancelled by other instanton amplitudes. Our manifestly 
supersymmetric scheme tells us that the most dominant contribution comes from Fig. 2 
including three insertions of pm so { t . The resulting correct Vhe is then estimated to be 



HE 



mis 



(_L)ia e -=fe Mp det(X u X d X h )\fy 2 (S) 2 A 
6 x 10 15 det{X u X d X h )\f X* 2 (GeV) 4 , 



(39) 




(40) 
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which is smaller than 10 _9 /^m 2 even when all Yukawa couplings are taken to be of order 
unity. In fact, in some cases one can draw an instanton graph without any insertion of 
P m soit- For instance, for the case with N h = 3 and n — I — k — 0, Fig. 3 provides an 
instanton amplitude which is consistent with all the internal symmetry selection rules but 
does not contain any pm so { t . Once again, such amplitude should be cancelled and in the 
manifestly supersymmetric scheme the dominant contribution is from Fig. 4 including two 
insertions of pm so f t . 

Let us finally comment a technical point about the number of loops in the instanton 
graph. When two graphs have the same value of N s , the dominant contribution comes 
from the one with a smaller number of loops. For instance, for the case with N h — 9 and 
(n, /, k) = (1, 1, 0), both (N 1: N 3 , N 4 , N 5 ) = (-3, 3, 0, 0) and (-3, -1,-1, -1) give the same 
value of Ns = 3/2. Fig. 5 shows the diagram of (—3,3,0,0) with the number of loops 
L = 11, while Fig. 6 is for (— 3, — 1, — 1, — 1) with the number of loops L = 13. In this 
example, Fig. 5 has two more scalar VEVs ((S)) than Fig. 6, but has two smaller number of 
loops. Usually a graph with more scalar VEVs has a fewer scalar loops because the number 
of scalar loops can be reduced by closing the scalar field line by inserting its VEV. Thus 
when there are two or more graphs with the same Ns, the graph with a more scalar VEVs 
has a fewer loops and thus gives a dominant contribution. 

To close all the quark zero modes in the instanton graph, one needs the insertion of 
either the bare quark mass or the Yukawa coupling involving a generic complex scalar field 
0. In our case, there is no bare quark mass and thus all the quark zero modes are closed by 
the Yukawa couplings. As a result, the instanton graph involves many 0-lines which should 
be closed again to yield a vacuum amplitude. There are many ways to close the 0-lines. 
One may close one 0-line without yielding any scalar loop by the insertion of (0), or close 
two 0-lines by inserting a complex mass of the form M|0 2 + h.c. which would yield one 
scalar loop. Inserting the couplings involving higher powers of to close the 0-lines results 
in more scalar loops. Obvioulsy the insertion of (0) in the instanton graph leads to the 
suppression by p(<f>). Also in supersymmetrc models, generically M? pa m so f t M^ where 
is a supersymmetric mass of the fermionic partner ip of 0, and thus the insertion of M| leads 
to a suppression by pm soft . Because of these suppressions, although economical in reducing 
the scalar loops, instanton graphs in which most of the scalar lines are closed either by (0) 
or M? do not give a dominant contribution to Vhe- This is the reason why in our cases 
the dominant contribution comes from the instanton graph involving relatively many scalar 
loops. 

4. CONCLUSION 

In this paper, we have studied the high energy axion potential Uhe induced by small QCD 
instantons in supersymmetric models in which QCD is not asymtotically free at high energy 
scales. In such models, Vhe may be larger than I0~ 9 f%m^ and then the axion solution to 
the strong CP problem can not be successfully implemented. To avoid the difficulty arising 
from the possible SUSY cancellation between different instanton graphs, we introduced a 
manifestly supersymmetric scheme to compute Uhe in generic SUSY models and applied it 
for a class of illustrative models. In our scheme, small instanton effects can be summarized 
by the effective superpotential (W e g) and also by the effective Kahler potential (K cS ). By 
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imposing the selection rules and a sensible limiting behavior together with the holomorphy, 
we showed that W e g does not allow a term which would contribute to the axion potential 
and thus Vhe always comes from K cS . As a simple consequence, Vhe is suppressed at 
least by two powers of pm so f t where p denotes the small instanton size and m soft is the soft 
SUSY breaking mass which is presumed to be of order the weak scale. In addition to this 
model-independent suppression, Vhe can be suppressed further by (i) small field VEVs (p4>i), 
(ii) small Yukawa couplings, and also (iii) loop factors. These model-dependent suppression 
factors are carefully analyzed for a class of illustrative models to estimate the possible largest 
value of Vhe and the results are summarized in Table 2. Our analysis suggests that the axion 
solution to the strong CP problem is stable against the effects of small instantons in a wide 
class of SUSY models even when QCD is not asymtotically free at high energy scales. 
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TABLES 



TABLE I. Quantum numbers of superfields and spurions 





U(1)a 


U(l)x 


U(l)x> 




Q 


1 








1 


u , a 


1 


— 1 





1 


i i r 

h, h c 


1 








1 







1 








s 








1 







l2+2N h 


-6 





6 




-2 











z h 


-2 





-(fc+l) 










-2 


-(n+1) 


2 


z 2 








-(1 + 3) 


2 


d 2 e 











-2 
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TABLE II. The largest possible size of Vhe an d S = Vhe/10 9 / 2 m 2 in each model. If S < 1, 
the axion solution to the strong CP problem is untouched by small instantons. 



N h 


n 


I 


k 


(N U N 3 ,N 4 ,N 5 ) 


N s 


L(loop) 


^he (GeV) 4 


5 


3 











(-3,0,0,0) 





8 


1.6 x 10~ 17 


4.1 x 10~ 5 


3 





1 





(-3,0,0,0) 





8 


1.8 x 10" 24 


4.5 x 10~ 12 


3 





1 


1 


(-3,3,0,0) 


3/2 


8 


1.8 x 10" 24 


4.5 x 10~ 12 


3 





2 





(-3,0,0,0) 





8 


4.1 x 10~ 32 


1.0 x 10~ 19 


3 





2 


1 


(-3,3,0,0) 


1 


8 


4.1 x 10~ 32 


1.0 x 10~ 19 


3 





2 


2 


(-3,1,-1,-1) 


4/3 


12 


2.7 x 10~ 28 


6.6 x 10~ 16 


3 


1 


1 





(-3,-3,0,0) 


3/2 


8 


2.0 x 10~ 63 


5.1 x 10~ 51 


3 


1 


1 


1 


(-3,0,0,0) 





11 


2.6 x 10~ 22 


6.6 x 10~ 10 


3 


1 


2 





(-3,-3,0,0) 


1 


8 


2.0 x IO- 47 


5.1 x 10~ 35 


3 


1 


2 


1 


(-3,0,0,0) 





11 


1.3 x 10~ 21 


3.3 x 10~ 9 


3 


1 


2 


2 


(-3,3,0,0) 


1 


11 


1.3 x 10~ 21 


3.3 x 10~ 9 


3 


2 


2 





(-3,-1,1,1) 


4/3 


12 


6.6 x 10" 75 


1.7 x 10~ 62 


3 


2 


2 


1 


(-3,-3,0,0) 


1 


11 


6.6 x 10^ 53 


1.7 x 10~ 40 


3 


2 


2 


2 


(-3,0,0,0) 





14 


4.3 x 10~ 27 


1.1 x 10" 14 


6 





1 





(-3,3,0,0) 


3/2 


8 


1.8 x 10" 24 


4.5 x 10" 12 


6 





2 





(-3,3,0,0) 


1 


8 


4.1 x 10~ 32 


1.0 x 10~ 19 


6 





2 


1 


(-3,4,-1,-1) 


7/3 


12 


2.7 x 10- 28 


6.6 x 10~ 16 


6 


1 


1 





(-3,0,0,0) 





11 


2.6 x 10~ 22 


6.6 x 10~ 10 


6 


1 


2 





(-3,0,0,0) 





11 


1.3 x 10~ 21 


3.3 x 10~ 9 


6 


1 


2 


1 


(-3,1,-1,-1) 


4/3 


15 


8.6 x 10~ 18 


2.2 x 10~ 5 


D 


o 
z 


o 
z 


n 
u 


{—0, —6, u, UJ 


1 
1 


1 1 


0.0 X 1U 


1.1 x iu 


6 


2 


2 


1 


(-3,-3,0,0) 


1 


14 


4.3 x 10~ 27 


1.1 x 10" 14 


9 





1 





(-3,2,-1,-1) 


2 


11 


5.9 x 10~ 14 


1.5 x 10- 1 


9 





2 





(-3,1,-1,-1) 


4/3 


12 


2.7 x 10~ 28 


6.6 x 10~ 16 


9 


1 


1 





(-3,3,0,0) 


3/2 


11 


2.6 x 10~ 22 


6.6 x 10~ 10 


9 


1 


2 





(-3,3,0,0) 


1 


11 


1.3 x 10~ 21 


3.3 x 10~ 9 


9 


2 


2 





(-3,0,0,0) 





14 


4.3 x 10~ 27 


1.1 x 10~ 14 
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FIGURES 



FIG. 1. Instanton graph which is consistent with all the internal symmetry selection rules but 
is required to be cancelled by SUSY for the case of = 9 and (n, I, k) = (0, 1, 0). Here the solid 
lines with and without waves denote the gluino and quark modes, respectively, while the dotted 
lines are the scalar fluctuations of Higgs, squarks, and S. The dark blobs represent the insertions 
of complex couplings or field VEVs which are explicitly written in the graph. 
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FIG. 3. Instanton graph which is required to be cancelled by SUSY for the case of Nh = 3 and 
(n,l,k) = (0,0,0). 




FIG. 4. Dominant instaton graph for the case of = 3 and (n,l,k) = (0,0,0). 
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FIG. 5. Instanton graph with the minimal values of N$ and the number of loops L for the case 
of N h = 9 and (n, I, k) = (1, 1, 0): N s = 3/2 and L = 11. 
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